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ABSTRACT 

We consider filters for the detection and extraction of compact sources on a background. We make a one¬ 
dimensional treatment (though a generalization to two or more dimensions is possible) assuming that the sources 
have a Gaussian profile whereas the background is modeled by an homogeneous and isotropic Gaussian random 
field, characterized by a scale-free power spectrum. Local peak detection is used after filtering. Then, a Bayesian 
Generalized Neyman-Pearson test is used to define the region of acceptance that includes not only the amplifi¬ 
cation but also the curvature of the sources and the a priori probability distribution function of the sources. We 
search for an optimal filter between a family of Matched-type filters (MTF) modifying the filtering scale such 
that it gives the maximum number of real detections once fixed the number density of spurious sources. We have 
performed numerical simulations to test theoretical ideas. 
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1. INTRODUCTION 

The detection of compact signals (sources) embedded in a background is a recurrent problem in many fields of 
science. Some common examples in Astronomy are the separation of individual stars in a crowded optical image, 
the identification of local features (lines) in noisy one-dimensional spectra or the detection of faint extragalactic 
objects in microwave frequencies. 

Regarding the detection of point sources on maps of the cosmic microwave background radiation (GMB), 
several techniques based on different linear filters have been proposed in the literature: the Mexican Hat Wavelet 
(MHW Cayon et al.,^ Vielva et al.^’^), the classic matched filter (MF, Tegmark and de Oliveira-Costa^), the 
Adaptive Top Hat Filter (ATHF, Chiang et al.®) and the scale-adaptive filter (SAF, Sanz et ah,® Herranz et 
al.^“®). A certain deal of controversy has appeared about which one, if any, of the previous filters is optimal for 
the detection of point sources in GMB data. 

In order to answer that question it is necessary to consider first a more fundamental issue, the concept of 
detection itself. The detection process can be posed as follows: given an observation, the problem is to decide 
whether or not a certain signal was present at the input of the receiver. The decision is not obvious since the 
observation is corrupted by a random process that we call ‘noise’ or ‘background’. 

Formally, the decision is performed by choosing between two complementary hypotheses: that the observed 
data is originated by the background alone (null hypothesis), and the hypothesis that the observation corresponds 
to a combination of the background and the signal. To decide, the detector should use the totality of the available 
information in terms of the probabilities of both hypotheses given the data. The decision device separates the 
space TZ of all possible observations in two disjoint subspaces, ^Z^: and TZ-, so that if an observation y G TZ- the 
null hypothesis is accepted, and if j/ G 7?., the null hypothesis is rejected, that is, a source is ‘detected’. Hence, 
we will call any generic decision device of this type a detector. 

Any detector can produce two kinds of errors: on the one hand, it can produce a false alarm or spurious 
detection when an observation in which no source was present is assigned to the subspace 7?,*. The probability 
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of this kind of error depends on the statistical properties of the background and the choice of the detector. On 
the other hand, a signal that is present in the observation can be missed by the detector (i. e. the observation 
is wrongly assigned to the subspace TZ-). This error is often referred as false dismissal. The probability of false 
dismissal depends on the statistical properties of the background, the choice of the detector and the properties 
of the signal (for example, its intensity). In general, it is not possible to decrease the incidence of both types of 
error at the same time: one of them can be reduced at expense of increasing the other. The goodness of a given 
detector must be established by taking into account the balance between these two types of error. 

The most simple example of detector, and one that has been exhaustively used in Astronomy, is thresholding. 
Thresholding considers that the space TZ of observations consists of all the possible values of the measured 
intensity f (in the case of an astronomical image) and subdivides this space into two simple regions TZ- = {^ G 
TZ: f < .J*} and 7^* = {^ G 7^ : ^ > ^h.}. The value is an arbitrarily chosen threshold that is often expressed as 
a number of times the standard deviation of the background, Thresholding works on the assumption 

that the probability of finding a value of f due to the background decreases as the value of f increases. In the 
case of a Gaussian background, this assumption has a very precise meaning and it allows us to straightforwardly 
control the probability of occurrence of spurious detections simply by setting a large enough threshold. However, 
this may lead to a very high probability of false dismissals. 

Unfortunately, in many cases the sources are very faint and this makes very difficult to detect them: a high 
threshold means that the number of detections will be very small. Here is where filtering enters in scene. The 
role of filtering is to transform the data in such a way that a detector can perform better than before filtering. 
For example, a filter can be designed in order to reduce the fluctuations of the background so that we can 
safely use lower detection thresholds and, hopefully, increase the number of detections without increasing the 
number of spurious detections. We remark that a filter is not a detector: the decision device we call ‘detector’ 
can be applied after the application of any imaginable filter, or even no filter at all, while the use of any filter 
without a posterior detection criterion means nothing. However, the two different steps in the process (filtering 
and detection) are not independent. In the thresholding example, the use of a filter that cancels most of the 
fluctuations in the background allows us to change the detection threshold from its original high value to a lower 
one. Given an adopted detector and a background, it is licit to ask which is the filter that creates the most 
favorables conditions in the filtered background for the detector to perform. In other words, the ‘optimality’ of 
a filter for detection depends on the type of detector chosen which, in turn, depends on the specific goal the 
observer has in mind: in certain cases the observer will accept a relatively large number of spurious detections 
in order to have a large number of true targets, whereas in other cases it could be more important to be certain 
that the detections are all of them reliable, and so on. 

For example, let us consider that we have chosen thresholding as our detection device. In that case, it has 
been shown that the optimal linear filter is the matched filter. It produces the maximum amplification of the 
signal with respect to the background fluctuations, so that the threshold for a given probability of spurious 
detections is minimum, allowing the thresholding detector to find more sources than would be detected if we 
filtered with any other filter. A sub-optimal approach, using only thresholding, is to select a priori a filter and 
adapt its scale in order to produce a maximum amplification in a given background. Such is the case of the 
Mexican Hat Wavelet at the optimum scale (MHO, Vielva et al.^’^) and the Adaptive Top Hat Filter (ATHF, 
Chiang et al.®). 

Thresholding has a number of advantages, among them the facts that it is straightforward, it has an obvious 
meaning in the case of Gaussian backgrounds, and it has been successfully used for many years in many fields of 
science. It, however, does not use all the available information contained in the data to perform decisions. The 
question is then: is it possible to devise a detector that uses additional information apart from mere intensities 
and that produces better results than thresholding? And, if so, which is the filter that optimizes the performance 
of such a detector? 

Let us focus on the case of one-dimensional data (such as stellar spectra, or time-ordered scannings of the sky 
in GMB experiments) and linear filters. Data in a one-dimensional array is entirely described by two quantities, 
namely the position in the array (corresponding to the spatial or temporal coordinate, for simplicity we will 
refer it as .spatial information) and the value (intensity) at each position. Thresholding uses only the intensity 
distribution to make the decision. Glearly, the inclusion of spatial information in a detector should be useful. 



For example, it could help to distinguish the sources from fluctuations in the background with similar scale but 
a different shape. A full description of this ‘spatial information’ should include the probability distribution of 
events (both due to background and sources) in space, with all its infinite moments. We will somewhat relax this 
demand of information assuming that the background is homogeneous and isotropic, and asking at each point 
for some information about the shape of the sources and the autocorrelation of the background (for example, 
the curvature of the peaks). 

In fact, even a simple filtering-and-thresholding scheme uses implicitly some degree of spatial information. 
Both the MHO and the ATHF adapt to the scale at which the contrast between sources and background produces 
the maximum amplification. The MF includes as well the information on the profile of the sources in order 
to amplify the structures whose shape correlates with the shape of the sources. SAF goes a step further in 
constraining additionally the scale of the filter. Moreover, in most cases the detection is performed not in all the 
points of the data but only in the peaks, that is, in those points where the curvature is positive. 

In a recent work, Barreiro et al.^*^ propose a detection criterion based on the Neyman-Pearson decision rule 
that uses the information of both the intensities ^ and the curvatures k of the peaks in a data set. In that work 
the performances of several filters (SAF, MF and MHW) is compared in terms of their reliability, defined as the 
ratio between the probability density of true detections over the probability density of spurious detections. They 
find that, on the basis of this quantity, the choice of the optimal filter depends on the statistical properties of the 
background. For the case of backgrounds that can be described with a power spectrum of the form P{q) oc q~'^, 
the SAF outperforms the other two filters for the case 1 < 7 < 1.6, whereas in the range 0 < 7 < 1 the MF is 
the most reliable. The MHW is the most reliable filter in this sense when 7 > 1.6. 

The reliability, defined in the previous sense, could not be a valid measure of performance of the Alters, 
since it favors a situation in which the number of the detections is very low in order to keep a ‘safe’ number of 
spurious detections. A different approach can be used in which the number density of spurious detections (or, 
alternatively, of true detections) is fixed for all the filters, and then the number density of true detections (or 
spurious detections) is compared for all the Alters. In this work, we first clearly define the goal of our experiment: 
fixed a certain number density of spurious detections, to obtain the maximum possible true detections from the 
data. Note that this goal is not universal: in other applications, an observer may desire to work on the basis of 
the reliability described in Barreiro et al.,^*^ or define its own requirements. Once fixed the goal, we will develop 
a detection criterion that corresponds to a Bayesian Generalized Neyman-Pearson test. Then, a particular 
filter from the family of matched-type filters (MTF) that optimizes the performance of the detector under the 
optimality conditions set by our goal will be obtained. Finally, the performance of the optimal MTF will be 
compared with the standard MF using simulations. 

The overview of this paper is as follows: In section 2, we introduce two useful quantities: number of maxima 
in a Gaussian background in the absence and presence of a local source. In section 3, we introduce the detection 
problem and define the region of acceptance. In section 4, we obtain different analytical and numerical results 
regarding point sources and scale-free background spectra and compare the performance of a new family of 
matched-type filters. In section 5, we describe the numerical simulations performed to test some theoretical 
aspects and give the main results. Finally, in section 6 , we summarize the main results and applications of this 
paper. 


2. BACKGROUND PEAKS AND COMPACT SOURCES 
2.1. The background 

Let us assume a ID background (e. g. one-dimensional scan on the celestial sphere or time ordered data set) 
represented by a Gaussian random field ^(x) with average value {S,{x)) = 0 and power spectrum P{q), q=\Q\'. 
(C(Q)C(*50) = P{<l)^D{q — q'), where ^(Q) is the Fourier transform of ^{x) and 6d is the ID Dirac distribution. 
The distribution of maxima was studied by Rice^^ in a pioneer article. The expected number density of maxima 
per intervals (x, x + dx), (v, v -|- dv) and (k, k -\- dn) is given by 
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where v S (— 00 , 00 ) and n G (0,oo) represent the normalized field and curvature, respectively, cr^ is the moment 
of order 2n associated to the field. 6c, dm are the coherence scale of the field and maxima, respectively. 


If the original field is linear-filtered with a circularly-symmetric filter R, b), dependent on 2 parameters 
{R defines a scaling whereas b defines a translation) 
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we define the filtered field as 

w{R,b) = j dx^{x)'^{x\R,b). (5) 

Then, the moment of order n of the linearly-filtered field is 
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being P{q) the power spectrum of the unfiltered field and tp{Rq) the Fourier transform of the circularly-symmetric 
linear filter. 


2.2. The presence of a local source 

Now, let us consider a Gaussian source (i. e. profile given by r(a;) = exp(—a;^/2i?^)) embedded in the previous 
background. Then, the expected number density of maxima per intervals {x, x + dx), [v, v P dv) and (k, k + dn), 
given a source of amplitude A in such spatial interval, is given by (Barreiro et al.^*^) 
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where v G (— 00 , 00 ) and k G (0,oo), Vg = A/ao is the normalized amplitude of the source and Kg = —AT'^la2 is 
the normalized curvature of the filtered source. The last expression can be obtained as 
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We consider that the filter is normalized such that the amplitude of the source is the same after linear filtering: 
f dx T(x)'i'(x; R, b) = 1. 


3. THE DETECTION PROBLEM 

We want to make a decision between filters based on detection. To make such a decision, we will focus on the 
following two fundamental quantities: a) the number of spurious sources which emerge after the filtering and 
detection processes and b) the number of real sources detected. As seen in the previous section, these quantities 
are properties of the Gaussian field and source that can be calculated through equations (1) and (7). As we will 
see, the previous properties are not only related to the signal-to-noise-ratio gained in the filtering process but 
depend on the filtered momenta to 4th-order (in the ID case), i. e. the amplification and the curvature of the 


source. 



3.1. The region of acceptance 

Let us consider a local peak in the ID data set characterized by the normalized amplitude and curvature (vg, Kg). 
Let Hq : n.d.f. nb{v, n) = n(v, k| 0) represents the null hypothesis, i. e. the local number density of background 
maxima, and Hi : n.d.f. n{h',K\h's) represents the alternative hypothesis, i. e. the local number density of 
maxima when there is a compact source with normalized amplitude and curvature {i>s,Kg = VgUg). Given the 
data {v, k), we can associate to any region k) two number densities nl and n* 
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n* = J p{h'g)di'g J n{h', K\h'g)di'dn. (10) 

Then, nl is the number density of spurious sources, i. e. due to the background, expected inside the region 
TZ^.{v,k), whereas n* is the number density of maxima expected in the same region of the {v,k) space in the 
presence of a local source. Hereinafter, we will call it the number density of detections. We will assume a 
Bayesian approach: at a concrete pixel we get the number of source detections weighting with the a priori 
probability p(vg). 7?,* is called the acceptance region. We remark that in order to get the true number of real 
source detections such a number must be multiplied by the probability to have a source in a pixel in the original 
data set. 

We will assume a Bayesian Generalized Neyman-Pearson decision rule using number densities instead of 
probabilities: the acceptance region 7^* giving the highest number density of detections n*, for a given number 
density of spurious nl, is the region (criterion for detection) 
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where L* is a constant. The proof follows the same approach as for the standard Neyman-Pearson test. Therefore, 
the decision rule is expressed by the likelihood ratio: if Tf > L* we decide that the signal is present, whereas if 
< L* we decide that the signal is absent. 

Once we have assumed the previous decision rule for detection, the region of acceptance 7?,* is given by 
L{v, k) > L* or equivalently by the sufficient linear detector 


7?.* : (p{v, k) > 
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where is a constant and p is given by 
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We remark that the assumed criterion for detection leads to a linear detector p (i. e. linear dependence on the 
threshold n and curvature k). 


3.2. Spurious sources and real detections 

Once obtained the region of acceptance 7^* in the previous subsection, one can calculate the number density of 
spurious sources and the number density of detections as given by equations (9) and (10). 
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Then, one can invert the equation for the number of spurious to get v?* = Pj Vs) that allows to rewrite the 

equation for the number of detections as n* = g{nl;0m, p,ys)- 


4. ANALYTICAL AND NUMERICAL RESULTS 


4.1. Point sources 


We will considerer as application the detection of compact sources characterized by a Gaussian proHle t{x) = 
exp(—a;^/2i?^), though the extension to other profiles will be considered in the future. Such a proHle is physically 
and astronomically interesting because represents the convolution of a point source (Dirac S distribution) with 
a Gaussian beam. 

4.2. The matched filter (ME) 

By introducing a circularly-symmetric Hlter, 4'(a:; B, b), we are going to express the conditions in order to obtain 
a matched filter for the detection of the source s(a:) at the origin taking into account the fact that the source 
is characterized by a single scale Ro- The following conditions are assumed: (1) {w{Ro,0)) = s(0) = A, i. e. 
w{Ro, 0) is an unbiased estimator of the amplitude of the source; (2) the variance of w{R, b) has a minimum at 
the scale i?o, i- e. it is an efficient estimator 
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This will be called matched filter as is usual in the literature. 
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For the case of a Gaussian proHle for the source and a scale-free power spectrum given by P{q) oc g the 
previous formula leads to the following matched filter 
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For the MF the parameters 0c, dm, P and the curvature of the source j/g are given by 
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We remark that the linear detector ip{v, k) is reduced to 
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i.e. the curvature does not affect the region of acceptance for the MF and the sufficient detector ip reduces to 
the plain thresholding detector. 

4.3. The family of Matched-type filters (MTF) 

Let us modify the filtering scale as aR to introduce a family of Matched-type filters. In particular, if we consider 
the white noise case (7 = 0 ), this family is given by 
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Therefore, this family allows one to filter at scales different from the one of the source. Obviously, for a = I the 
usual MF is recovered. 

The parameters characterizing the background and source are in this case modified in the following way 
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4.4. Uniform distribution of point sources 

In this case, 

P(i^s) = 

Vc 

This allows one to obtain 
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4.5. Theoretical results 

Hereinafter, we shall consider the case of white noise (7 = 0) for the background and a uniform distribution for 
the point sources in an interval ly € [ 0 , i^c]- has been chosen to produce a threshold of 2 for the hltered held 
using the standard MF. The corresponding thresholds for the MMF can be easily obtained using equation (23). 
We would like to point out that this reflects a situation where we want to detect very weak sources, which were 
below the Ict level in the original unhltered map. 

The region of acceptance 7?.* is dehned by the sufficient linear detector ip. Let us consider equation (14), 
giving the number density of spurious, then one can invert to get (p*(n^). In Fig. 1 we show some results for 
the cases R=3, nl = 0.05 and i?=4, nl = 0.03. The two panels on the left correspond to the hrst case and the 
ones on the right to the second one. The top figures show the expected number density of detections versus the 
a parameter of the hlter. It is interesting to remark that such a number decreases with a. In fact, the largest 
number of detections are obtained for a ~ 0.3 which corresponds approximately to filter at the pixel scale. The 
relative ratio to the standard MF dehned as r = (n*(a)/n*(a = 1) — 1) x 100 is shown in the bottom panels. For 
the hrst case, the ratio takes values up to ~ 10 per cent whereas for the second it goes up to ~ 5 per cent. These 
results clearly indicate that, under certain conditions, the standard MF can be improved by simply modifying 
the scale of the hlter. We remark the importance of the curvature dehning the acceptance regions for the MTF. 


5. NUMERICAL SIMULATIONS: RESULTS 

5.1. The simulations 

In order to test the previous ideas, we have simulated a set of one-dimensional images containing a Gaussian 
background characterized by a white noise power spectrum (7 = 0 ) and point sources distributed in intensity 
following an uniform distribution such as in eq. (24). Since we are interested in detecting very faint sources, we 
set the upper-limit threshold cut Vc = 2 with respect to the MF. Specihcally, we simulated a set of images with 
4096 pixels each, with white noise dispersion unity (background image) and then added a source with a Gaussian 
prohle of FWHM=7 pixels, i.e. i? = 3, at the central pixel of each image (background-|-source image). The size 
of the image is such that the addition of the source does not modify the previous dispersion in a significant way. 
Then, each image was filtered with the MTF given by equation (21). 

For consistency, all the relevant quantities needed for the detection were estimated directly from the images. 
The observables are the amplitude and curvature of the maxima as well as the moments cr". From these observed 
quantities, it is possible to calculate the parameters p and j/s of the filtered images and, therefore, the value of 
the sufficient detector p (see equation 13) for each peak. The value yj* that dehnes the acceptance region can be 



R =3, rij,* 0.05 


R = 4, nb* = 0.03 






Figure 1. The expected number density of detections (top panels) and the relative ratio to the standard MF (bottom 
panels) are given for the MTF for different values of the a parameter. Two cases are considered: R = 3, nl — 0.05 (left 
figures) and R = A, nl = 0.03 (right figures). 

obtained from the images as well following a very simple procedure. For any background maxima in the image, 
the detector ip is calculated. Then all the values of p are sorted. Once we fix the number of spurious sources 
nl, the value of (/?* is given by the p above which the number of maxima found corresponds to nl. Once defined 
the acceptance region, we can proceed to apply the detector to any maxima where the existence of a source is 
suspected. In our simulations, we looked at the central pixel of the images were the sources were introduced and 
applied the detection criterion every time a maximum was found in that pixel. 

5.2. Results 

The results of the simulations are shown in figure 2 for the case i? = 3, nj=0.05, where we represent the number 
of detections n* vs. the a parameter of the filter. The dotted line represents the result obtained from the 
simulations whereas the solid line is the theoretical result. We have done five simulations for each a value to 
estimate the error bars. In order to have 10000 images with a maximum of the background in the central pixel, 
for each of the simulations we generate between ~ 50000 and ~ 150000 realizations (depending on a). The two 
numbers reflect the fact that for lower a the same nl is achieved with a smaller number of realizations. We 
remark the agreement between the two results and also that ~ 10 per cent more sources are found filtering at 
the pixel scale {a ~ 0.3) than with the standard MF (a = 1). Therefore, we have shown that in conditions of 
white noise and for weak sources it is possible to outperform the number of detections given by the MF using 
the detection criterion presented here and an appropriate MTF. This clearly suggests the use of other optimal 
filters under certain conditions, different from the MF, from the practical point of view, for the detection of weak 
sources. 

This methodology is optimal for weak sources embedded in white noise. Regarding bright sources, the 
amplification plays a major role as compared to curvature entering the acceptance region. Therefore, the use of 
the MF is well justified in this case. 

Regarding the behavior of the detector with respect to the scale R, we found that if R is very low the sources 
become more and more point-like and they are more easily mistaken with the fluctuations of the background. 








Figure 2. Number density of detections for the case R = 3, nl — 0.05 for different values of the a parameter. The solid 
line represents the expected theoretical value whereas the dotted line is obtained through the numerical simulations. 


Therefore, in these cases the curvature does not help to distinguish among them and the detector tends again 
to the thresholding case which favors the use of the MF. 

6. CONCLUSIONS 

The detection of compact sources on a background is a relevant problem for many fields, in particular for Astron¬ 
omy. Several detection techniques based on the use of linear filters thresholding-based detectors are standard. 
Here, we have considered an approach to the problem of detector design based on a Bayesian generalization of 
the Neyman-Pearson rule that includes a priori information of the source distribution and the number densities 
of maxima (background and background plus source) to dehne the acceptance region. 

Our approach based on maxima includes both the amplitude and the curvature. Therefore, the chances of 
detection do not depend only on the amplification of the sources produced by the hltering but also on the hltered 
momenta up to the fourth order. This determines in a strong way the designing of the linear hlters that are used 
to help the detection. We have applied our technique to a family of matched-type hlters (MTF) by modifying 
the scale of the standard matched hlter. We considered the very interesting case of white noise to represent the 
background. As an example, we have considered a uniform distribution of sources in the interval 0 < v < 2 
in the hltered held, i.e. weak sources. We have shown that the curvature plays an important role dehning the 
acceptance region and we have proven that the number of detections in the case of a hlter with a scale similar 
to the pixel size beats the number of detections in the case of the standard MF. This result has been tested with 
numerical simulations for a uniform distribution and white noise. 

The ideas presented in this paper can be generalized: application to other prohles (e. g. multiquadrics, 
exponential) and non-Gaussian backgrounds is physically and astronomically interesting. The extension to 
include several images (multi-frequency) is relevant. The generalization to two-dimensional data sets (plane and 
spherical maps) and nD images is also very interesting. Finally the application of our method to other helds is 
without any doubt. We are currently doing research in some of these topics. 





6.1. Acknowledgments 

The authors thank Enrique Martmez-Gonzalez and Patricio Vielva for useful discussions. RBB thanks the 

Ministerio de Ciencia y Tecnologia (MCYT) and the Universidad de Cantabria for a Ramon y Cajal contract. 

DH acknowledges support from the European Community’s Human Potential Programme under contract HPRN- 

CT-2000-00124, CMBNET. MLC thanks the MCYT for a predoctoral FPI fellowship. We acknowledge partial 

support from the Spanish MCYT projects ESP2001-4542-PE and ESP2002-04141-C03-01 and from the EU 

Research Training Network ‘Cosmic Microwave Background in Europe for Theory and Data Analysis’ 

REFERENCES 

1. Cayon, L., Sanz, J. L., Barreiro, R. B., Martmez-Gonzalez, E., Vielva, P., Toffolatti, L., Silk, J., Diego, 
J. M. & Argiieso, F., ’’Isotropic wavelets: a powerful tool to extract point sources from cosmic microwave 
background maps”, MNRAS 315, pp. 757-761, 2000. 

2. Vielva, P., Martmez-Gonzalez, E., Cayon, L., Diego, J. M., Sanz, J. L. & Toffolatti, L., ’’Predicted Planck 
extragalactic point-source catalogue”, MNRAS 326, pp. 181-191, 2001. 

3. Vielva, P., Barreiro, R. B., Hobson, M.P., Martmez-Gonzalez, E., Lasenby, A., Sanz, J. L. & Toffolatti, L., 
’’Combining maximum-entropy and the Mexican hat wavelet to reconstruct the microwave sky”, MNRAS 
328, pp. 1-16, 2001. 

4. Tegmark, M. & de Oliveira-Costa, A., ’’Removing Point Sources from Cosmic Microwave Background Maps”, 
ApJL 500, pp. 83-86, 1998. 

5. Chiang, L.Y., Jprgensen, H.E., Naselsky I.P., Naselsky P.D., Novikov I.D. & Christensen, P.R., ”An adaptive 
filter for the construction of the Planck Compact Source Catalogue”, MNRAS 335, pp. 1054-1060, 2002. 

6 . Sanz, J. L., Herranz, D. & Martmez-Gonzalez, E., ’’Optimal Detection of Sources on a Homogeneous and 
Isotropic Background”, ApJ 552, pp. 484-492, 2001. 

7. Herranz, D., Sanz, J. L., Barreiro, R. B. & Martmez-Gonzalez, E., ’’Scale-adaptive Filters for the Detec¬ 
tion/Separation of Compact Sources” ApJ 580, pp. 610-625, 2002. 

8 . Herranz, D., Sanz, J. L., Hobson, M., Barreiro, R. B., Diego, J. M., Martmez-Gonzalez, E. Lasenby, A.N., 
’’Filtering techniques for the detection of Sunyaev-Zel’dovich clusters in multifrequency maps” MNRAS 336, 
pp. 1057-1068, 2002. 

9. Herranz, D., Gallegos, J., Sanz, J. L. & Martmez-Gonzalez, E., ’’Point source detection and extraction from 
simulated Planck time-ordered data using optimal adaptive Hlters”, MNRAS 334, pp. 533-541, 2002. 

10. Barreiro, R. B., Sanz, J. L., Herranz, D. & Martmez-Gonzalez, E., ’’Comparing filters for the detection of 
point sources”, MNRAS 342, pp. 119-133, 2003. 

11 . Rice, S. O., ’’Mathematical Analysis of Random Noise” in ’’Selected papers on Noise and Stochastic Pro¬ 
cesses”, N. Wax, ed., Dover Publ. Inc. (N. Y.) 1954, pp. 133-294. 



